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■^H NMR spin-lattice relaxation and line-shape analyses are performed to study the temperature- 
dependent dynamics of water in the hydration shells of myoglobin, elastin, and collagen. The results 
show that the dynamical behaviors of the hydration waters are similar for these proteins when using 
comparable hydration levels of ^ = 0.25 — 0.43. Since water dynamics is characterized by strongly 
nonexponential correlation functions, we use a Cole-Cole spectral density for spin-lattice relaxation 
analysis, leading to correlation times, which are in nice agreement with results for the main dielectric 
relaxation process observed for various proteins in the literature. The temperature dependence can 
roughly be described by an Arrhenius law, with the possibility of a weak crossover in the vicinity 
of 220 K. Near ambient temperatures, the results substantially depend on the exact shape of the 
spectral density so that deviations from an Arrhenius behavior cannot be excluded in the high- 
temperature regime. However, for the studied proteins, the data give no evidence for the existence 
of a sharp fragile-to-strong transition reported for lysozyme at about 220 K. Line-shape analysis 
reveals that the mechanism for the rotational motion of hydration waters changes in the vicinity of 
220 K. For myoglobin, we observe an isotropic motion at high temperatures and an anisotropic large- 
amplitude motion at low temperatures. Both mechanisms coexist in the vicinity of 220 K. ^^C CP 
MAS spectra show that hydration results in enhanced elastin dynamics at ambient temperatures, 
where the enhancement varies among different amino acids. Upon cooling, the enhanced mobility 
decreases. Comparison of ^If and ^'^C NMR data reveals that the observed protein dynamics is 
slower than the water dynamics. 



I. INTRODUCTION 

The biological function of many proteins depends on 
the existence of a hydration shell [l| . Despite consider- 
able progress in recent years, a complete understanding 
of the underlying coupling of protein and water motions 
is still lacking. In order to obtain fundamental insights, 
it proved useful to study the temperature dependence of 
protein and water dynamics, in particular through the 
protein dynamic transition at about 200 K ji] , which 
was found to be related to a freezing of biological func- 
tions [J]. The protein dynamics were proposed to be 
"slaved" [1| or " plasticized" by water dynamics, but 
the nature of the interplay of protein and water dynamics 
is still a subject of controversial debate. 

Traditionally, protein dynamics were characterized on 
the basis of temperature-dependent atomic mean square 
displacements obtained from neutron scattering (NS) or 
Mossbauer experiments 2]. Then, a crossover from lin- 
ear behavior to non-linear behavior indicates an onset 
of non-harmonic motion upon heating, which was taken 
as evidence for the protein dynamic transition. More re- 
cently, combined NS and dielectric spectroscopy (DS) ap- 
proaches were used to determine the temperature depen- 
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dence of correlation times r for protein dynamics p, |7| ■ 
Since a discontinuity was not observed it was argued that 
a true dynamic transition does not exist, but the phe- 
nomena result from the protein dynamics entering the 
experimental time windows when the temperature is in- 
creased [1, 

Various experimental methods were employed to inves- 
tigate the temperature-dependent dynamics of protein 
hydration waters. NS studies found a sharp crossover 
from a non- Arrhenius (fragile) behavior to an Arrhenius 
(strong) behavior upon cooling and attributed it to a 
fragile-to-strong transition (FST) at ca. 220 K ^'l . In- 
terestingly, the postulated FST occurs in the vicinity of 
the protein dynamical transition and a relation of both 
phenomena was postulated j!,, J,0], in harmony with re- 
sults from molecular dynamics (MD) simulations 
However, other studies, including DS work, did not ob- 
serve a sharp transition, but a weaker crossover at lower 
temperatures or no crossover at all [H, [13, [3j IHj [3 ■ 
These authors concluded that water dynamics, at least 
at sufficiently low temperatures, is not dominated by the 
structural (a) relaxation, but a possibly largely univer- 
sal [IB| secondary (/3) relaxation. For example, water 
dynamics was related to the Johari-Goldstcin (JG) /? re- 
laxation p^ . 

For a fundamental understanding of the interplay of 
protein and water dynamics, it is necessary, first, to as- 
sign the observed relaxation processes to the components 
and, second, to identify them as a or /3 relaxations, 
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requiring determination of the mechanisms for the re- 
spective motion. NMR experiments and MD simulations 
proved very useful tools for this purpose. In recent stud- 
ies fisf . we exploited the capabilities of these tech- 
niques to determine the temperature-dependent mecha- 
nisms for the dynamics of water in the hydration shells of 
elastin and collagen, two main components of the connec- 
tive tissue. It was found that the mechanism for the wa- 
ter motion continuously changes from an isotropic motion 
to an anisotropic motion upon cooling. Moreover, the lat- 
ter low-temperature process was shown to involve large- 
angle jumps. Such anisotropic large-angle jump motion 
is at variance with the mechanisms for the a and JG-/3 
relaxations of molecular supercooled liquids, suggesting 
the existence of a water-specific (3 process. However, ex- 
act origin and possible universality of this /3 process are 
open questions. 

Here, we perform temperature-dependent ^H and ^^C 
NMR experiments on hydrated proteins to investigate 
the water and the protein dynamics, respectively. On 
the one hand, we analyze ^H NMR spin-lattice relaxation 
(SLR) and ^H NMR hue shape (LS) to study the respec- 
tive water dynamics in the hydration shells of myoglobin, 
elastin, and collagen. In this way, we intend to address 
the question of universality. Moreover, for myoglobin, we 
perform a detailed comparison of the present results with 
literature data [l^ to ascertain whether NMR and DS 
yield a coherent picture of water dynamics. On the other 
hand, protein dynamics is investigated in ^^C NMR using 
cross polarization (CP) and magic-angle spinning (MAS) 
techniques. Specifically, we record ^^C CP MAS spectra 
of elastin to monitor site-resolved changes of the protein 
dynamics upon hydration and cooling, respectively. 



II. MATERIALS AND METHODS 
A. Materials and experiments 

We study collagen, elastin, and myoglobin hydrated to 
different levels h (g water per 1 g protein) using D2O (^H 
NMR) or H2O (i^C NMR). The samples are referred to 
as C (collagen), E (elastin) or M (myoglobin) followed by 
the hydration level in percent. Details about the colla- 
gen and elastin samples can be found in Ref. M35 
was prepared from sperm-whale myoglobin (Sigma) and 
kindly provided by W. Doster. 

^H NMR SLR and LS measurements were performed 
at Larmor frequencies luo of 2tt x 76.8 MHz for collagen 
and elastin samples and 2tt x 46.1 MHz for M35. The for- 
mer setup is described in Ref. [l^]. In the latter setup, a 
home-built probe was placed in a Konti CryoVac cryostat 
operating with liquid nitrogen and controlled by a TIC 
304 MA CryoVac temperature controller. The tempera- 
ture stability was better than 0.5 K. A 90° pulse length 
of 2.1 /is was used, ensuring adequate excitation of the 
broad spectrum. The solid-echo pulse sequence, 90° - tp 
- 90°, preceded by a saturation-pulse sequence, was em- 



ployed to acquire ^H NMR spectra. The ^H NMR SLR 
was measured utilizing the saturation-recovery technique 
in combination with a solid echo for detection. The ^^C 
CP MAS NMR spectra were recorded with a Bruker DSX 
400 spectrometer working at ojq — 27t x 100.6 MHz and 
a 4-mm triple-resonance Bruker MAS probe. For CP, a 
^H 90° pulse of 3.5 fis length was followed by a contact 
time of 1 ms. Two-pulse phase modulation (TPPM) de- 
coupling [2l[ with a field strength of 71 kHz was applied. 
The spinning speed was 8 kHz. The ^^C chemical shifts 
were referenced with respect to tetramethylsilane. A flow 
of nitrogen gas controlled by a Bruker VT 3000 heat- 
ing unit was used to adjust sample temperature. Since 
set and actual temperature differ substantially in MAS 
experiments, deviations were removed within an uncer- 
tainty of ±1 K by means of temperature calibration with 
2"7pb NMR of lead nitrate [H. 



B. Theoretical background 

In solid-state ^H NMR, the dominant interaction is the 
first-order quadrupolar interaction, which describes the 
coupling between the electric quadrupole moment of the 
deuteron and the electric field gradient (EFG) at the nu- 
clear site. The corresponding quadrupolar contribution 
to the NMR frequency is given by [23j : 

wq(6',0) = ±- [3cos^{e) - 1 -r;sin2(6')cos(20)]. (1) 

Here, 9 and are the polar coordinates of the external 
magnetic flux density Bq in the principal axis system 
of the EFG tensor, S and rj are the anisotropy parameter 
and the asymmetry parameter, respectively, and the "±" 
sign corresponds to the two allowed transitions of the 
1 = 1 nucleus. In our case of D2O, the asymmetry of the 
EFG tensor is negligibly small (77 « 0.1), leading to 

uq{9) = ±^-[3cos^9)-1], (2) 

where 9 is the angle between Bq and the direction of the 
0-D bond. Thus, the orientation of the water molecule 
determines the resonance frequency, rendering ^H NMR 
a precious tool to study water rotational dynamics. 

If the correlation time of water reorientation ex- 
ceeds the inverse anisotropy parameter, t ^ 1/S, the 
anisotropy of the quadrupolar interaction and the powder 
average will result in broad ^H NMR spectra for disor- 
dered materials ^25] . When the temperature is increased 
and, hence, molecular dynamics becomes faster, these 
static spectra, which are called Pake spectra, collapse as 
soon as t becomes comparable to 1/6, i.e., for molecular 
dynamics in the microseconds regime. Finally, motion- 
ally averaged spectra are observed in the limit of fast mo- 
tion, T ^ 1/6. In particular, fast isotropic reorientation 
results in narrow Lorentzian lines. When a broad distri- 
bution of correlation times G(logr) exists so that there 
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are both slow (r » 1/5) and fast (r <C 1/^ ) molecules at 
a given temperature, the LS will be comprised of a Pake 
spectrum and a motionally narrowed line [11] . Such spec- 
tra have been called two-phase spectra. 

For isotropic motions, the SLR time Ti is given by 

funs 



1 
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[J(c^o) + 4J(2c^o)] 



(3) 



where J(w) is the spectral density of the molecular re- 
orientation. Hence, SLR analysis provides access to the 
time scale of molecular rotational motion, e.g., Ti is a 
minimum when r « I/wq, i.e., for molecular dynamics in 
the nanoseconds range. 

Quantitative SLR analysis requires knowledge about 
the shape of the corresponding rotational autocorrelation 
function F2{t). For a Debye process, F2{t) = exp(— f/r). 
Then, the spectral density can be written as 



Jbpp{uj) 



1 + (ujTy 



(4) 



and Eq. ([3]) takes the form derived by Bloembergen, Pur- 
cell, and Pound (BPP) [27||. For disordered materials, ex- 
istence of distributions G(log t) usually leads to a nonex- 
ponentiality of F2{t). In the literature [l^jISll, different 
spectral densities were used to consider such distribu- 
tions. For various hydrated proteins, DS studies found 
that the loss peak, which was ascribed to water dynamics, 
is symmetric and well described by the emp irical Cole- 
Cole (CC) distribution function 0,[ilIli,[ll[il- There- 
fore, we use the CC spectral density [2g] 



sin(f acc)(wTcc)"'^ 



w[l + (wrcc)2"" -1-2 cos(f acc)(wacc)" 



(5) 



in our SLR analysis. Here, Tcc is the position of the 
maximum of G(logT) and Ucc describes the width of this 
symmetric distribution. Near ambient temperatures, in- 
formation about the shape of the spectral density from 
DS is lacking. In this temperature range, we also employ 
the Cole-Davidson (CD) spectral density 



sm[Pcd arctan(tJTerf)] 
u;[l + (a;Ted)2]/3c./2 



(6) 



for ^H SLR analysis. The asymmetric CD distribution 
is characterized by a mean correlation time Ted — TcdPcd 

M. 



III. RESULTS 

A. NMR SLR 

Figure [T] shows a typical buildup of the normalized 
equilibrium magnetization M{t) after saturation for M35. 
We see a two-step shape, which is a consequence of 
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FIG. 1: Buildup of the normalized equilibrium magnetization 
after saturation for M35 at 231K. The line is a fit with Eq. (O 



deuteron-proton exchange between water and myoglobin 
molecules [1^ : exchangeable hydrogens of the protein are 
partially replaced by deuterons so that the ■^H NMR sig- 
nals are comprised of the respective contributions from 
water deuterons and myoglobin deuterons. The obser- 
vation of distinct steps means that the deuteron-proton 
exchange is slow on the time scale of the buildup and, 
hence, we can separately study the dynamical behaviors 
of the components in ^H NMR. In previous work on sam- 
ples with higher hydration levels of /i > 0.6 [l3|, we ob- 
served a third step at long times, which was attributed 
to crystalline water. For M35, we find no evidence for 
the existence of such third step and, thus, of crystalline 
water. 

Considering the two-step signature of the buildup, we 
fit the data to 



M{t) = 1 



+ {I - aw)e 



I3p' 



(7) 



Here, Ti^ and Tip are the SLR times of water and pro- 
tein, respectively, and and (1 — a^) describe the 
heights of the steps. To take into account that ^H NMR 
SLR relaxation can be nonexponential for disordered ma- 
terials [1^, we introduce the stretching parameters (3^ 
and f3p and calculate the mean SLR times as 



{Tiw,ip) 



Ti 



w,lp 



1 



(8) 



where F is the Gamma function. The indices w and p 
denote water and protein, respectively. 

The results of this analysis enable unambiguous assign- 
ment of the steps to water and protein, see Fig. [2l The 
SLR time of one component is short and and exhibits 
a pronounced minimum, whereas that of the other com- 
ponent is long and increases upon cooling in the whole 
studied temperature range. While the latter behavior in- 
dicates a slow (r ^ l/i^o) component and, hence, the 
long-time step can be assigned to the protein, the former 
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FIG. 2: a) Arrhenius plot of the mean SLR times (Ti^) 
and (Tip) of water and protein deuterons in M35, respec- 
tively. (Ti) is the mean SLR time in the temperature range 
T < 155 K where one-step magnetization buildup is ob- 
served. The experimental results are compared with SLR 
times {Tiw)r)S, which were calculated using the distributions 
G(logr) from DS on M45 [l^]- Inset: stretching parameter 
characterizing the nonexponentiality of the SLR step from 
water deuterons. b) Temperature dependence of the relative 
heights a™ of the water signals for M35, E25, E43, and C25. 



behavior is characteristic of a highly mobile component 
so that the short-time step results from the water. This 
assignment is confirmed by the relative heights of the re- 
laxation steps at ambient temperatures. On the basis 
of the amino acid composition of myoglobin [30| . it is 
expected that water deuterons substantially outnumber 
myoglobin deuterons for M35. Therefore, the higher first 
and lower second steps result from water and protein, re- 
spectively. The pronounced temperature dependence of 
ttw is an NMR effect, as will be discussed below. 

Having assigned the steps, the SLR of water and myo- 
globin can be discussed in more detail. For water, we 
find a minimum of (Tiu,) at about 240 K, indicating water 
dynamics with r « I/wq- Important information about 
a possible existence of a distribution G'(logT) is offered 
by the value at the minimum, {Ti^)m- It is straightfor- 
ward to calculate {Tiw)m for a Debye process, see Eq. 
gD. Using ujQ = 2n X 46.1MHz and S = 2tt x 166kHz, 
as determined from the low-temperature spectrum in 



Fig. [4l a minimum value of 1.4 ms is expected, which is 
about a factor of 2 smaller than the experimental result 
{Tiw)m = 3.6ms. This discrepancy means that a Debye 
process does not describe the water dynamics, but a dis- 
tribution G'(logT) exists [2j]. For myoglobin, (Tip) shows 
a weak increase with decreasing temperature, implying 
T ^ l/'i^o- Below 160 K, equilibrium magnetization is 
build up in a single step and, thus, water and myoglobin 
exhibit a common SLR time of a few seconds. This means 
that an exchange of magnetization occurs between the 
components on this time scale. Possible candidates for 
exchange processes are spin diffusion, i.e., a transfer of 
magnetization due to flip-flop processes of the spins, and 
deuteron exchange. 

Next, we turn to the stretching parameters j3w and 
/3p, describing deviations from exponential SLR. For wa- 
ter, we find exponential SLR above 200 K, while cooling 
results in a decrease of the stretching parameter [3^ at 
lower temperatures, indicating increasingly nonexponen- 
tial SLR, see Fig. [2ja). In our case, nonexponential SLR 
is expected since the existence of a distribution of corre- 
lation times should lead to a distribution of SLR times, 
V{Ti). Therefore, our findings imply that, above 200 K, 
water dynamics averages over any distribution V{Ti)^ 
i.e., each water molecule samples a representative set of 
local environments on the time scale of Ti^ , leading to ex- 
ponential SLR. When the temperature is decreased, the 
time scale of water dynamics and, thus, of exchange pro- 
cesses within the distribution G(logr) approaches that 
of SLR so that the averaging process becomes incom- 
plete and the existence of a distribution V{Ti) is revealed, 
i.e., SLR becomes nonexponential. Such temperature de- 
pendence of the stretching parameter is known from ^H 
NMR SLR studies on the a process of supercooled liquids 
[23 |. For myoglobin, the stretching parameter /3p = 0.6 
is independent of temperature within the accuracy of the 
measurement. Thus, the protein dynamics are hetero- 
geneous and the corresponding distribution V{Ti) is not 
averaged due to exchange processes within the distribu- 
tion G'(logT). 

Before we continue our SLR analysis, it is instructive 
to determine whether NMR and DS studies yield con- 
sistent results for the dynamics of myoglobin hydration 
water. For this purpose, we directly use the CC distribu- 
tions obtained from a DS study on M45 [l^ to calculate 
the SLR times according to Eq. (O. The resulting val- 
ues {Titu)ds are included in Fig. [2fa). We see that the 
calculated data nicely agree with the NMR results in the 
common temperature range below the (Ti^) minimum. 
Thus, NMR and DS probe the same dynamical process, 
which is appropriately described by a CC distribution, at 
least below 220 K. 

For further SLR analysis, we first assume that a CC 
distribution describes water dynamics at all studied tem- 
peratures. In DS [l^, an essentially temperature inde- 
pendent width parameter ttcc ~ 0.48 was found for M45. 
In ^H NMR, the width parameter is determined by the 
value of (Ti^,)m. We obtain — 0.45 at 240 K, in 
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FIG. 3: Correlation times of protein hydration waters. The 
values of Tcc, Ted, and Taw from the present SLR analyses 
for water in the hydration shells of myoglobin (M35), elastin 
(E25), and collagen (C25) are compared with literature data 
for myoglobin and lysozym (L) hydration waters. Specifically, 
correlation times from DS on M45 19], NS on M35 [S^l, DS 
on L37 [?[, and NS on L30 ^9] are shown. For comparison, 
correlation times of protein dynamics from NS on L40 Q are 
included as line. 



harmony with the result from DS. Assuming that acc 
does not depend on temperature, we can extract Tcc from 
(Tik,). In Fig. [21 we compare the temperature depen- 
dence of Tcc for hydrated myoglobin, elastin and collagen. 
The values for the two latter proteins are obtained when 
we use the CC spectral density {acc — 0.50) to analyze 
the SLR times presented in previous work JJl]. We see 
that the hydration waters of all studied proteins exhibit 
similar correlation times, which roughly follow an Arrhe- 
nius law. Closer inspection of the data implies that there 
may be a weak change of the temperature dependence 
at about 220 K. Consistently, separate fits of the myo- 
globin data in the ranges T > 220 K and T < 220 K yield 
activation energies of Ea — 0.66 eV and Ea — 0.55 eV, 
respectively. 

Since information from DS about the shape of the spec- 
tral density is lacking for myoglobin above 220 K, we de- 
termine to which extent SLR analysis depends on the 
choice of J(w) in this temperature range. Using Jcd{'^), 
a width parameter (3cd = 0.19 results from the value of 
{Tiw)m- Assuming a temperature independent width of 
the CD distribution, we determine the mean correlation 
times fed from Tiyj . The results are included in Fig. [3l 
We see that the CC and CD spectral densities yield com- 
parable correlation times in the range 230 — 270 K, i.e., 
in the vicinity of the (Ti^,) minimum, while the findings 
differ substantially near ambient temperatures. In par- 
ticular, non-Arrhenius behavior is observed for Ted- 

Figure [2l^b) shows the temperature-dependent relative 
height of the water signal in M35 together with data 
for E25, E43 and C25. For the studied samples, 
hibits a comparable temperature dependence with a min- 



imum at about 205 K. This minimum can be explained 
by the use of a solid echo for signal detection in our ex- 
periments, see Sec. Ill Al It indicates that water dynamics 
during the dephasing and rephasing periods tp — 20 fis 
of the solid-echo sequence interferes with echo formation 
[3l| . This conclusion is confirmed by the results of our 
SLR analysis, which shows that water dynamics is indeed 
characterized by r of the order of a few microseconds at 
the minimum position of 205 K, see Fig. [3l By contrast, 
the signal is not reduced for much faster and slower dy- 
namics when ujQ is time- averaged and time- invariant, re- 
spectively. Detailed random-walk simulations for differ- 
ent motional models show that the solid-echo arnplitude 
is a minimum for a correlation time of about 3 /xs (isl . Is^ . 
Hence, the coincidence of the minima of shows that 
the hydration waters of myoglobin, elastin and collagen 
exhibit a correlation time Taw = 3/is at 205 K. In Fig. 
[31 we see that Taw is comparable to Tcc, indicating that 
both methods yield consistent results. 



B. NMR line-shape analysis 

In order to gain information about the mechanisms for 
water dynamics, we analyze the temperature-dependent 
NMR LS. Results for M35 are presented in Fig. [2 
Fully relaxed (FR) and partially relaxed (PR) spectra are 
shown in the left and right columns, respectively. When 
measuring FR spectra, we wait for complete recovery of 
magnetization after saturation, while the waiting time 
between the saturation and solid-echo sequences was set 
to {Tiw) when acquiring PR spectra. Thus, both myo- 
globin and water contribute to the FR spectra, while the 
water component is singled out in the PR spectra. Com- 
parison of the FR and PR spectra shows that, at all stud- 
ied temperatures, myoglobin adds a broad Pake spectrum 
to the LS, indicating that the deuterated regions of the 
protein show rotational motion that is slow (r ^ 1/6) 
or has small angular amplitude (< 15°). Moreover, we 
see that the relative contribution of myoglobin to the FR 
spectra is particularly prominent in the vicinity of 205 K, 
where the contribution of water is reduced because of dy- 
namics during the solid-echo sequence, see Fig. [2Jb). 

Focusing on the PR spectra, we see that, below 200 K, 
the spectral contribution of water is comprised of a broad 
Pake spectrum and a central feature, the prominent in- 
tensity of which extends ca. from —20 to 20 kHz. With 
decreasing temperature, the intensity of the latter com- 
ponent slowly decreases, while its shape does not change. 
Unfortunately, it is not possible to single out the water 
signal in PR spectra below 186 K because the difference 
between Tiw and Tip diminishes. However, it is evident 
from the FR spectra at low temperatures, that the central 
feature continues to loose intensity upon further cooling 
until it disappears between 100 and 150 K. The obser- 
vation that, below 200 K, the central part of the spec- 
trum changes intensity, but not shape shows that it is 
a fast-motion limit spectrum, i.e., it results from water 
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FIG. 4; Solid-echo spectra of M35 at various temperatures. 
Left column: fully relaxed spectra comprising signals from 
myoglobin and water deuterons. Right column; partially re- 
laxed spectra to single out the signal from water deuterons, 
see text for details. 



molecules exhibiting t <^l/5. The shape of this compo- 
nent clearly deviates from a Lorentzian, indicating that 
the underlying motion is anisotropic. The coexistence 
of spectra typical of the limits of fast and slow motion 
implies that a broad distribution G(log r) governs the 
anisotropic water reorientation at T < 200 K. Then, the 
slow (r ^ 1/^) water molecules from the distribution 
give rise to the Pake spectrum, while the fast [t <^l/6) 
ones yield the central feature between —20 and 20 kHz. 
When temperature is varied and, hence, G(logT) shifts, 
the relative intensity of both spectral components con- 
tinuously changes, while the respective LS does not [26j . 

Above 205 K, a narrow Lorentzian line emerges, indi- 
cating that a fraction of water molecules start exhibit- 
ing fast (r <Si 1/^) isotropic motion, or at least fast 
motion with cubic symmetry. Interestingly, both fast- 
motion limit spectra, the Lorentzian and the central fea- 
ture between —20 and 20 kHz coexist in the tempera- 
ture range 208 — 218 K. Hence, there are water molecules 
that exhibit fast isotropic motion while others show fast 
anisotropic motion, but no fast isotropic motion. This 



means that two different mechanisms exist for the ro- 
tational motion of hydration water. Upon heating, the 
intensity of the Lorentzian rapidly increases until the PR 
spectra are well described by a single Lorentzian above 
230 K. Hence, more and more water molecules exhibit 
fast isotropic motion and, eventually, all molecules are 
involved in such dynamics in the high-temperature range. 

The shape of the central feature provides important 
insights into the mechanism for water rotational motion 
below 200 K. The main intensity of this LS component 
extends from —20 to 20 kHz and minor intensity is found 
up to about ±40 kHz. Thus, as compared to the static 
case, the width of the spectrum is reduced by more than 
a factor of 4. Such pronounced reduction of the spectral 
width cannot result from small-amplitude motion, but 
it rather requires large-amplitude motion. Furthermore, 
two-site jumps lead at most to a motional narrowing by a 
factor of 2 ^23]. Hence, the small width of the central fea- 
ture provides strong evidence against the possibility that 
the water molecules perform exact tt flips about their 
6*2 symmetry axes. For fast uniaxial rotation or for fast 
jumps between N sites {N > 3) with Cjy symmetry, the 
motionally averaged quadrupolar coupling tensor is char- 
acterized by an anisotropy parameter ^ = 0. Moreover, 
since 77 w in the static case, the anisotropy parameter 
5 of the averaged tensor is given by [2^ 

^ = i(3cos2/3-l). (9) 

Here, f3 is the angle between the rotation axis and the 
possible orientations of the 0~D bonds. Then, our case of 
\S\ « 27r X 40 kHz translates into /3 « 45° or /3 « 65°. The 
ambiguity results because the symmetry of the LS ham- 
pers determination of the sign of the asymmetry parame- 
ter in ^H NMR. Thus, within these models, the maximum 
angular displacement of the 0-D bonds, which corre- 
sponds to the full opening angle 2/3 of the cone, amounts 
to 90° or 130°. 

However, the observed central feature cannot be de- 
scribed by a single Pake spectrum, but it is comprised 
of a distribution of spectral shapes, corresponding to a 
distribution of geometries of water motion, which may 
be expected for a disordered material. For example, a 
fraction of the motionally averaged tensors may exhibit 
asymmetry parameters 77 > 0. Therefore, the present 
results are at variance with exact tt flips, but not with 
strongly distorted tt flips. Also, LS analysis allows us 
to rule out exact tetrahedral jumps, which would lead 
to a narrow Lorentzian in the limit of fast motion, but 
strongly distorted tetrahedral jumps cannot be excluded. 
Altogether, the shape of the central feature does not en- 
able determination of a unique motional model, but it 
nevertheless indicates that, below 200 K, myoglobin hy- 
dration water performs anisotropic reorientation, which 
covers large angular displacements, but deviates from ex- 
act two-site or tetrahedral jumps. 

Finally, we study whether the ^H NMR LS differs for 
the hydration waters of various proteins. FR spectra of 
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E25 are shown in Fig. [3 We see that the hmit of slow 
motion, t 3> l/i5, is found somewhat below 145 K. As the 
temperature is increased, a central, broad feature gains 
intensity and, finally, the LS is dominated by a narrow 
Lorentzian. Thus, the results for E25 resemble that for 
M35, see Fig.^l However, the shape of the central feature 
is different for both proteins. In particular, for E25 and 
E43 [3 at T 210 - 220 K, the data give no evidence 
for a coexistence of spectral components indicative of fast 
isotropic and fast anisotropic motions, respectively. The 
origin of this difference is not yet understood. It may 
indicate that the mechanisms for low-temperature ro- 
tational motion differ between the hydration waters of 
elastin and myoglobin. Alternatively, a crossover from 
isotropic to anisotropic motion may occur at larger val- 
ues of r for elastin than for myoglobin so that the limit 
of fast motion cannot be reached for the anisotropic dy- 
namics of elastin hydration waters. 



FIG. 7: "C CP MAS NMR spectra of E60 at various tem- 
peratures. 



C. "C CP MAS NMR spectra 

To investigate the dependence of elastin dynamics on 
hydration and temperature, we use ^'^C NMR. Figures [5] 
and [7] show temperature-dependent ^^C CP MAS spec- 
tra of EO and E60, respectively. The data were divided 
by the number of used accumulations and corrected for 
the Curie factor, enabling straightforward comparison. 
For EO, the results are independent of temperature in 
the studied range 186 - 298 K, indicating that substan- 
tial elastin dynamics is absent on the experimental time 
scale. By contrast, a strong temperature dependence of 
the spectra is evident for E60. At 298 K, the spectral 
intensity is significantly reduced and the relative inten- 
sities of the various peaks differ from that found for EO. 
In ^•^C CP MAS, such reduction of the peak intensities 
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results when molecular dynamics in the sub- milliseconds 
regime reduce the strength of the ^H-^^C dipolar cou- 
plings, which are required for effective polarization trans- 
fer. Therefore, the LS of E60 shows an enhancement of 
elastin dynamics due to hydration near room tempera- 
ture, consistent with findings in previous works |34l. Iss!] . 
The strength of this effect differs for different types of 
carbons. It is prominent for Gly Cq. and Pro Ca,p,-y, im- 
plying that the dynamics of Gly and Pro amino acids are 
particularly sensitive to hydration. When the tempera- 
ture is decreased, the LS of E60 continuously approaches 
the temperature-independent LS of EO and, eventually, 
the spectra of both samples are very similar at T < 243 K. 
Thus, the enhanced elastin mobility of the wet sample 
diminishes. Consistently, the results of the present and 
previous [17.] NMR studies on hydrated myoglobin, 
elastin and collagen revealed an increase of (Tip) upon 
cooling, see Fig. [21 implying a continuous reduction of 
protein mobility. 



IV. DISCUSSION 

We have used and -'^^C NMR to study temperature- 
dependent water and protein dynamics, respectively. 
NMR LS and SLR analyses have revealed that, despite 
minor quantitative differences, water shows comparable 
rotational dynamics in the hydration shells of collagen, 
elastin, and myoglobin. For all studied materials, water 
reorientation cannot be described as a Debye process, but 
broad distributions of correlation times G'(logT) exist. 

To consider the existence of such distributions, SLR 
analysis has mainly been performed using the CC spec- 
tral density. This choice has been motivated by the obser- 
vation of symmetric dielectric loss peaks for various hy- 
drated proteins 0, [13, [3 • The mean correlation times 
Tec and width parameters dec obtained from such SLR 
analyses are in good agreement for the hydration waters 
of collagen, elastin, and myoglobin. The temperature 
dependence of Tcc is roughly described by an Arrhenius 
law, although there is some evidence for a weak crossover 
at about 220 K. The width parameters acc ~ 0.48 corre- 
spond to full widths at half maximum of about 3.3 orders 
of magnitude. Near ambient temperatures, DS does not 
provide information about the shape of the spectral den- 
sity. Therefore, we also employed the CD spectral density 
for SLR analysis. While the correlation times from Jcd 
and Jcc are comparable in the range 230 — 270 K, they 
differ at higher temperatures. In particular, application 
of the CD spectral density leads to deviations from Ar- 
rhenius behavior in the high-temperature range. 

LS analysis has revealed that the mechanism for rota- 
tional motion of protein hydration waters changes in the 
vicinity of 220 K, where the temperature-dependent cor- 
relation times may exhibit a weak crossover. The nature 
of these changes differ somewhat in the hydration shells 
of diverse proteins. For myoglobin hydration water, we 
have observed an isotropic motion at high temperatures 



and an anisotropic motion at low temperatures. Interest- 
ingly, both mechanisms coexist in the vicinity of 220 K, 
suggesting that two different dynamical processes exist. 
Moreover, LS analysis has shown that the anisotropic 
motion has large amplitude. Specifically, it involves an- 
gular displacements on the order of the tetrahedral an- 
gle. Our findings exclude that water performs exact two- 
site or tetrahedral jumps. However, distorted two-site 
jumps, including inexact tt flips, or distorted tetrahe- 
dral jumps may be consistent with the data, provided 
there is a substantial distribution of geometries of the 
motion. Although LS analysis does not enable deter- 
mination of a unique model of low-temperature water 
reorientation, it allows us to draw some important con- 
clusions. Because of its anisotropy, the observed motion 
is not related to the a process of myoglobin hydration 
water. Moreover, the large-amplitude of the reorienta- 
tion is strong evidence against the conjecture that 
the low-temperature process can be identified with the 
JG-/3 process, which results from small-amplitude mo- 
tion for molecular glass formers [38| . Therefore, we pro- 
pose that a water-specific (3 process is observed, consis- 
tent with the conclusions in our previous experimental 
and computational studies [13, [3 ■ Specifically, in MD 
simulation work, we found that the mechanism for wa- 
ter motion changes from small-angle (diffusive) motion to 
large-angle (jump) motion [18] . In future work, wc plan 
to use two-dimensional NMR in time and frequency 
domains to further narrow down the possible mechanisms 
for the rotational dynamics of protein hydration waters 
at low temperatures. 

It is instructive to compare the present results with 
literature data for the dynamics of protein hydration wa- 
ters. First, we emphasize that, although there may be 
a weak crossover in the temperature-dependent correla- 
tion times in the vicinity of 220 K, the NMR data of all 
studied hydrated proteins are at variance with an exis- 
tence of a sharp FST reported in NS work for the hy- 
dration water of lysozyme at such temperatures [§] , see 
Fig. [21 There, it was found that the translational re- 
laxation time of hydration water follows an Arrhenius 
law with a small activation energy of Ea — 0.14eV below 
220 K, consistent with temperature-dependent diffusivity 
measurements (36j . A dynamical process with such char- 
acteristics is not observed in our present and previous 
[l3] NMR studies on hydration waters. For example, 
if there was a FST, the low-tcmperature process would 
be the a process and, hence, fast isotropic motion with 
T <Sil/5 would exist down to about 140 K, see Fig.[3l In 
this case, a narrow Lorentzian is expected as NMR 
LS, at variance with the observations at such tempera- 
tures, see Figs. [Hand [3 

Below 220 K, NMR and DS correlation times nicely 
agree, see Fig. [31 indicating that the same dynamical 
process is probed. For both techniques, the depicted cor- 
relation times describe the peak position of G(logT) and, 
hence, direct comparison should be possible. However, 
we have to consider that DS and NMR probe the rota- 
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tional correlation functions of the first and second Leg- 
endre polynomials, respectively. The time constants of 
both these correlation functions coincide when the reori- 
entation involves large jump angles, while they differ by 
a factor of 3 for small jump angles, i.e., for rotational 
diffusion 2^. Hence, the observed agreement of NMR 
and DS data may suggest that hydration water shows 
large-angle jumps below 220 K. For the width of the CC 
distributions, comparable values acc ~ 0.48 are obtained 
in NMR and DS Thus, both methods yield consis- 
tent distributions G(logr). In particular, for myoglobin, 
the values of (Ti™) measured in NMR agree with that 
calculated on the basis of the distributions G'{logT) de- 
termined in DS work [l9l |. However, the correlation times 
of the present SLR analysis are about two orders of mag- 
nitude shorter than that determined in our previous ^H 
NMR stimulated-echo study . This apparent discrep- 
ancy can be rationalized when we consider that these 
NMR techniques measure different types of mean corre- 
lation times. While stimulated-echo experiments yield 
the mean correlation time (r) , which is the time integral 
of the rotational correlation function, the peak position 
Tec of G(logr) results from the present SLR analysis. For 
a logarithmic Gaussian distribution of correlation times 
characterized by a width parameter a, both time con- 
stants are related by [13] 

log.o (^) - ^ (10) 
\tccJ 2 logio e 

when exchange within G(logr) can be neglected. Us- 
ing a — 1.43 from the present SLR analysis, we expect 
logio ((■'")/'^cc) = 2.35, consistent with the experimentally 
observed difference of both mean values. 

Above 220 K, DS 0, [H, 13 and NS 0, H studies 
reported deviations from an Arrhenius law for various 
hydrated proteins and a Vogel-Fulcher-Tammann (VFT) 
equation was used to interpolate the data. In this tem- 
perature range, the present SLR analysis has suffered 
from lacking knowledge about the shape of the spectral 
density. While VFT behavior results from use of the CD 
spectral density, Arrhenius behavior is obtained from the 
CC spectral density. In any case, we have found that, al- 
though a distribution of correlation times exists, SLR 
is exponential above 220 K, indicating that the water 
molecules sample a representative set of local environ- 
ments on the time scale of Ti^,, i.e., exchange processes 
within G(logr) result in short lifetimes of the dynamical 
heterogeneities. Such behavior is known from ^H NMR 
SLR studies on the a process of supercooled liquids ^23] . 

Altogether, our results may suggest the following pic- 
ture of the temperature-dependent dynamics of protein 
hydration waters. Above about 220 K, water exhibits 
isotropic rotational motion, which is governed by pro- 
nounced dynamical heterogeneities with short lifetimes. 
Hence, the dynamics resembles the a process of su- 
percooled liquids, although the vicinity of protein sur- 
faces may lead to some differences, e.g., in the shape of 
G(logr) 0, [H, [il- Below about 220 K, water exhibits 



anisotropic large-amplitude rotational motion, implying 
existence of a water-specific [3 process, which shows sim- 
ilar features for the hydration waters of various materi- 
als, including polymers and inorganic nanoporous matri- 
ces jli, [l5j] . The crossover between both mechanisms is 
smooth, suggesting a direct relation of both processes. 
For example, the [3 process may be a spatially restricted 
precursor of the a process. Alternatively, both processes 
may require breaking of the same number of hydrogen 
bonds, leading to comparable temperature dependencies 
and distributions G(logT). For example, interstitial de- 
fect diffusion cannot be excluded as origin of the low- 
temperature process [13] • Because of the large angu- 
lar displacements, the water-specific (3 process destroys 
a large fraction of orientational correlation. Thus, only a 
small fraction is left for the a process, if existent [12], so 
that it is difficult to extract information about the struc- 
tural relaxation from NMR and DS studies below 220 K, 
say. Furthermore, it is an open question whether the ob- 
served change in the motional mechanism is related to 
an ordering of the hydrogen-bond network upon cooling, 
e.g., as a consequence of a proposed transition from a 
high-density liquid to a low-density liquid [39j . 



^•^C CP MAS spectra have shown that hydration re- 
sults in enhanced elastin dynamics at ambient temper- 
atures. This enhancement varies among different amino 
acids and it is particularly prominent for Gly and Pro. 
Upon cooling, the enhanced elastin mobility decreases 
and vanishes at about 243 K. The finding that LS changes 
cease at higher temperatures in ^^C NMR than in ^H 
NMR indicates that protein dynamics is slower than wa- 
ter dynamics. Combined application of DS and NS on hy- 
drated lysozyme enabled identification of two relaxation 
processes of the protein 0. While the slower protein 
process exhibits much longer correlation times, the time 
constants of the faster process coincide with that of the 
water dynamics, see Fig. [3] Therefore, the faster process 
was ascribed to a water-coupled structural relaxation of 
lysozyme. The present finding that protein dynamics is 
slower than water dynamics implies that the faster pro- 
tein process is not probed in the present ^'^C CP MAS 
spectra, suggesting that this process cannot be identified 
with the structural relaxation, although possible differ- 
ences between elastin and lysozyme need to be studied 
in future work. Rather, ^'^C CP MAS spectra reveal a 
freezing of elastin dynamics faster than milliseconds at 
about 243 K. At the present, the origin of the observed 
protein process is an open question. Moreover, it is not 
yet clear whether this process involves those degrees of 
freedom that freeze in during the protein dynamic tran- 
sition. Our finding that, at all studied temperatures, 
the proteins contribute a Pake spectrum to the ^H NMR 
LS shows that essentially all sub-milliseconds dynamics 
of deuterated protein regions have small amplitudes, in 
harmony with results in the literature [s^ . [3a ] , 
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